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1.0 


INTRODUCTION 


In  the  ABEL  device,  the  laser  beam  must  pass  through  an 
interface  between  two  flowing  gases  of  different  indices  of 
refraction,  the  main  laser  gas  and  nitrogen  in  the  end  flow,  as 
illustrated  in  Figure  1.  Inhomogeneous  turbulent  mixing  in  the 
interface  causes  phase  aberration  and  degradation  of  the  laser 
beam. 


If  the  amount  of  phase  aberration  \j,'^  is  small,  the  beam 
degradation  is  given  by: 


(1) 


where  I  and  IQ  are  the  degraded  and  nondegraded  beam  intensi¬ 
ties  respectively.  A  theoretical  expression  of  the  beam  degrada¬ 
tion  due  to  turbulence  has  been  derived; (2) 


(2) 


where  X  is  the  laser  wavelength,  fin  is  the  refractive  index 
difference  across  the  mixing  layer;  a  is  an  average  scalar 
fluctuation  intensity,  L  is  the  mixing  layer  thickness  and  l  is 
the  turbulent  macroscale  related  to  the  refractive  index  fluctua- 
t  ion. 


Experimental  studies  on  beam  degradation  across  a  turbulent 
mixing  layer  were  carried  out  in  a  small-scale  flow  test  in  which 
the  mixing  layer  was  observed  in  a  20  cm  high  optical  chimney. 


1.  Born,  M.  and  Wolf,  E. ,  Principles  of  Optics,  Perraaqon  Press, 
1975,  pp.  463-465. 

2.  Vu,  B.T.  ,  Sutton,  G.  ,  Theophanis,  G.  and  Limpaecher,  P.., 
"Laser  Beam  Degradation  Through  Optically  Turbulent  Mixing 
Layers,"  AIAA  Paper  80-1414,  13th  AIAA  Fluid  and  Plasma 
Dynamics  Conference,  Snowmass,  Colorado,  1980 
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The  four  specific  objectives  of  the  test  were: 

1)  to  obtain  experimental  verification  of  the  theoretical 
expression  for  beam  degradation  due  to  turbulence  (E)q.  2), 

2)  to  study  beam  degradation  characteristics  and  mixing 
layer  growth  under  velocity  sheared  conditions, 

3)  to  test  a  method  for  reducing  the  beam  degradation  of  a 
single  layer  by  reducing  the  mixing  layer  growth,  and 

4)  to  study  the  mixing  layer  stability  if  one  stream  is 
exposed  to  a  region  of  stagnant  fluid,  as  in  the  case  of 
the  end  flow  shown  in  Figure  1. 

Four  series  of  experiments  were  carried  out  in  connection 
with  the  four  above  objectives. 


in 
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2.0  TEST  APPARATUS 


2.1  EXPERIMENTAL  APPARATUS 

Experiments  at  room  temperature  were  carried  out  for  a  mixing 
layer  between  a  stream  of  nitrogen  and  a  stream  of  a  mixture  of 
helium  and  nitrogen.  Mach-Zehnder  interferometer  data  together 
with  temperature  and  velocity  data  were  obtained. 

The  flow  setup  is  shown  in  Figure  2  and  is  described  in 
detail  elsewhere.^)  The  two  streams  of  gases  divided  by  a 
splitter  plate  in  a  plenum  chamber,  passed  through  a  flow  plate 
and  were  allowed  to  mix  in  an  optical  chimney  7.11  cm  long, 

5.08  cm  wide  and  20.0  cm  high.  Two  types  of  flow  plates  were  used 
in  the  experiments  and  are  shown  in  Figures  3a  and  3b.  The  flow 
plate  shown  in  Figure  3a  is  identical  to  that  of  the  ABEL  laser; 
the  modified  form  shown  in  Figure  3b,  had  a  small-mesh  screen 
positioned  over  the  orifice  plate  whose  dimensions  are  shown  in 
Figure  3c.  Proper  construction  of  the  modified  flow  plate  is 
important  and  is  described  in  more  detail  in  Section  4.3.  The 
plenum  chamber  pressures  were  measured  by  two  pressure  gauges  with 
60  psig  maximum  pressure  and  0.2  psi  subdivisions.  The  pressures 
were  constant  within  0.5  psi  for  each  experimental  run.  The 
velocities  of  the  gases  in  the  optical  chimney  were  measured  with 
a  hot-wire  anemometer  and  were  also  cross  checked  against  predic¬ 
tions  from  the  plenum  pressures  as  discussed  in  Section  2.2. 
Temperatures  of  both  streams  were  measured  in  the  chimney.  They 
were  constant  for  all  gases  used  and  were  about  3°F  below  room 
temperature . 

2.2  CALIBRATION 

Gas  velocity  measurements  were  made  with  a  hot-wire  anemo¬ 
meter  (Thermal  Science,  Inc  1050}  and  a  standard  hot-film  sensor 
( TS 1/1210-20  Sensor,  0.002"  d  ia .  ,  0.06"  long ,  frequency  response 
250  kHz).  The  hot-film  sensors  were  calibrated  with  nitrogen 
flow,  helium  flow  and  five  mixtures  of  nitrogen  and  helium. 

Typical  calibration  curves  are  shown  in  Figure  4. 


3. 


"Cold  Flow  Electric  Laser  Development  {U)“,  ABEL  I 
Technical  Report  for  Period  1  Jan.  1977  -  30  Sept. 
MICDM  Contract  DAAK40 -7  5-C-i  272,  AERL  Doc.  79-320, 
1979.  CONFIDENTIAL 
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The  velocity  of  each  stream  measured  by  trie  hot-film  sensor 
averaged  across  the  flow  plate  was  correlated  with  the  plenum 
pressure.  Figure  5  shows  the  calibrated  data  for  different  gas 
mixtures  at  different  plenum  pressures.  The  solid  lines  in  the 
figure  represent  theoretical  velocities  calculated  from  the  plenum 
pressures.  The  theoretical  calculations  are  based  on  subsonic  and 
choked  flow  through  the  orifice  plate. 

For  subsonic  flow,  the  velocity  in  the  chimney  is 


V  = 


-  P) 
P 


RT 

M 


(3) 


where 


Cd  is  the  discharge  coefficient 

Aq/A  is  the  open-area  ratio  of  the  flow  plate 

PQ  is  the  plenum  pressure 

P  is  the  static  pressure  of  the  gas  in  the  optical 

chimney 

R  is  the  universal  gas  constant  (8314  J/kg  mole) 

M  is  the  molecular  mass  (kg/kg  mole) 

T  is  the  gas  temperature  ("K) 

If  the  plenum  pressure  is  larger  than  the  critical  pressure 


P  > 
o 


(P) 


(4) 


where  is  the  ratio  of  the  specific  heats  Cp/C v,  then  the 
sonic  choked  flow  relation  applies 


V  = 


P 

o 

P 


(5) 


where 

T0  is  the  total  temperature  of  the  gas  (°K) 

T  is  the  static  temperature  of  the  gas  in  the  optical 

c  h imney 


17 


AVCO  EVERETT 


^ZAVCO  EVERETT 


K  is  ( Y  +  1 )  /  ( >  -  1 ) 

To  simplify  calculations ,  we  take  Ac/A  as  the  open-area  ratio 
for  the  orifice  plate  alone  which  is  6.7%.  Flow  plate  calibration 
is  obtained  with  hot-wire  anemometer  measurements  (Fig.  5) ,  and 
the  nominal  discharge  coeffients  are  found  to  be: 

Cd  =  0.69  +  3.5%, 
for  the  screened  flow  plate, 
and  =  0.79  +  5%, 

for  the  unscreened  flow  plate. 
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3.0  TEST  DESCRIPTION 


3.1  PROCEDURE 

The  two  principal  test  measurements  were  the  beam  degradation 
measurement  and  the  mixing  layer  growth  rate.  These  two  measure¬ 
ments  must  be  related  to  the  gas  velocities  and  the  refractive 
index  difference  (or  gas  mixture  composition)  across  the  mixing 
layer. 

The  velocities  were  either  measured  by  hot-wire  anemometry  or 
calculated  using  Eq.  (3)  and  Eq.  (5)  with  the  measured  discharge 
coefficients.  For  hot-wire  measurements,  the  errors  were  about 
2%,  while  velocities  calculated  from  plenum  pressures  had 
uncertainies  of  about  5%. 

The  refractive  index  difference  was  measured  directly  from 
the  inte rf erograms ,  with  the  probe  beam  oriented  parallel  to  the 
mixing  layer.  The  fringe  shift  across  the  mixing  layer  was 
obtained  by  comparing  a  tare  inter ferogram  without  flow  and  a 
disturbed  interf erogr am  with  flow.  The  number  of  shifted  fringes 
A J  is  -elated  to  the  difference  in  refractive  index  An  as 

A  J  W  (6) 

where  x  is  the  laser  wave  length  (5260  ft)  and  w  is  the  width  of 
the  optical  chimney  in  the  optical  axis.  With  w  =  5.08  cm  for  our 
experiment , 

An  =  AJ  x  1.04  x  10~4 

If  the  mixing  layer  is  between  a  stream  of  nitrogen  and  a  stream 
of  a  mixture  of  nitrogen  and  helium,  both  at  temperature  T  (°K) 
and  1  atm.  ,  then  the  mixture  composition  and  the  refractive  index 
difference  are  related  by: 


An 


<V 


273 


(7) 


where  £>N2and  3^  are  the  G1  adstone- Da le  constants  for  nitrogen 
and  the  mixture  respectively.  The  constant  BM  is  given  by 


d 


M 


+  X 


He 


(8) 
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where  and  B ^  ace  respectively  the  molecular  fraction  and  the 
Gladstone-Dale  constant  of  species  "i".  Figure  6  shows  An  and  the 
normalized  velocity  of  the  N2/He  mixture  under  sonic  conditions 
as  a  function  of  composition. 

3.2  DATA  ANALYSIS 

The  inter ferogr am  data  reduction  was  done  by  the  following 
methods.  First,  the  mixing  layer  growth  was  studied  by  orienting 
the  probe  beam  parallel  to  the  layer.  If  the  profile  of  the 
average  refractive  index  n  in  the  mixing  layer  is  given  by  the 
error  function,  we  can  define  the  mixing  layer  thickness  L  as  the 
maximum  slope  thickness. (2)  Figure  7  illustrates  this  method  and 
shows  that  at  the  edqe  of  this  thickness,  n  differs  from  the  value 
of  n  in  the  free  stream  by  0.105  An.  From  the  disturbed  interfero- 
gram,  we  can  thus  determine  the  edge  of  maximum  slope  thickness  of 
the  mixing  layer.  This  method  of  determining  the  thickness  is 
more  accurate  and  objective  than  the  "visual"  determination  of  the 
mixing  layer  thickness  base  on  shadowgraphs  such  as  used  by  Brown 
and  Roskho.(^)  The  visual  method  usually  gives  an  overestimation 
for  a  mixing  layer  with  larger  \n  as  compared  to  a  mixing  layer 
with  small  An. 

Second,  the  value  of  beam  degradation  across  the  layer  is 
determined  from  interf erograms  taken  with  the  probe  beam  oriented 
perpendicular  to  the  mixing  layer.  It  is  convenient  to  deal  with 
the  square  root  of  the  beam  degradation  or  the  RMS  of  the  phase 
aberration,  yip '  2  ,  since  they  are  linearly  proportional  to  An  and 
the  thickness,  L,  of  the  mixing  layer.  The  RMS  phase  aberration  is 
obtained  by  the  following  method.  A  portion  of  the  disturbed 
interf erogram  is  first  digitized.  The  linear  dimensions  of  the 
analyzed  area  of  the  inter ferogr am  are  chosen  to  be  at  least  ten 
times  the  turbulent  macroscale.  With  the  fringes  numbered  consec¬ 
utively,  the  digital  data  are  the  x-y  coordinates  of  the 
intersection  of  the  fringes  with  the  reference  lines  drawn  nearly 
perpendicular  to  the  fringes,  the  distances  between  consecutive 
reference  lines  are  chosen  to  be  comparable  to  the  fringe 
spacings.  Computer  analysis  of  the  digital  data  were  done  with  a 
computer  code  developed  by  the  Optical  Sciences  Center  (OSC)  at  the 
University  of  Arizona.  A  detailed  description  of  the  OSC  code  is 
given  in  the  Appendix  A.  The  program  takes  the  (x,y)  array  points 
and  fits  them  to  a  polynomial  of  the  form 

f  <  x , y)  =  Ax  +  By  +  C  +  D(x2  +  y2)  (3) 


4.  Brown,  G.L.  and  Roshko,  A.,  "On  Density  Effects  and  Large 

Structure  in  Turbulent  Mixing  Layers,"  J.F.M.  Vol.  64,  1974, 
pp.  775-816. 


22 


.  AVCO  EVERETT 


}  {  erf  [o  ( x ) y]  -I  )  +  I  (see  ref.  2 ) 


Mixing  Layer  Thickness 
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using  the  least-square  method.  The  first  three  terms  represent  the 
tilt;  while  the  last  term  gives  the  defocus.  The  aberration  given 
by  f(x,y)  from  a  plane  is  deterministic  and  does  not  contain  any 
perturbation  due  to  turbulence.  The  difference  between  the  fit  and 
the  stored  array  of  data  yields  the  phase  errors,  i.e., 

OPDi  =  ORDi  -  f(xi,yi)  (1 

where  OPD^  is  the  optical  path  difference  at  point  i  in  the 
matrix,  and  ORD  is  the  fringe  order  at  point  i.  The  RMS  phase 
aberration  due  to  turbulence  is  then  computed  from  the  RMS  value  of 
OPD .  Using  this  analysis  we  can  thus  effectively  compute  the 
random  phase  aberration  due  to  turbulence  es  separated  from  any 
deterministic  aberration. 

As  mentioned  earlier,  within  each  sample  area  of  the 
i nte rf erograms ,  digital  samplings  were  made  at  intervals  comparable 
to  the  fringe  spacings,  which  were  about  2  mm  in  our  experiments. 
Since  the  turbulent  macroscale  varies  with  streamwise  distance,  the 
sampling  interval  can  be  larger  or  smaller  than  the  turbulence 
scale  size.  However,  in  calculating  the  RMS  phase  aberration,  the 
ratio  of  the  sampling  interval  to  the  turbulence  scale  size  has  no 
effect  on  the  accuracy  of  the  remits.  This  is  so  since  the  RMS 
pi  ^se  aberration  calculated  here  can  be  considered  as  obtained  by 
ensemble  averaging.  In  section  4.0  it  is  shown  that  the 

values  of  a  related  to  the  RMS  phase  aberration  obtained  in  these 
experiments,  throughout  the  shear  layers,  are  in  agreement  with 
values  obtained  in  all  previous  turbulent  measurements  in  shear 
flows.  The  agreement  confirms  the  adequacy  of  the  sampling 
technique  used  here. 

The  largest  uncertainty  in  the  RMS  phase  aberration  results  is 
due  to  reading  of  the  interferometric  fringes,  which  were  of  finite 
thickness.  The  thickness.,  as  observed  in  the  interf erogr am,  is 
about  0.3  of  the  fringe  spacings.  If  the  reading  errors  are  random 
about  the  center  of  the  fringe  an  estimated  error  in  phase  aberra¬ 
tion  is  given  as 


*  -  +  - *  —  =  +  0.10  fringe  spacing. 

NT 

Scatter  in  actual  readings,  as  presented  in  the  next  section, 
is,  however,  only  about  one  half  of  this  estimate,  since  care  was 
exercised  in  locating  the  center  of  each  fringe. 


5.  Bend at,  J.S.  and  Piersol,  A.G.,  Random  Data:  Analysis  and 
Measurement  Procedures,  J.  Wiley  &  Sons,  p.  68,  87,  288, 
(1971)  . 
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4.0 


TEST  DATA  AND  RESULTS 


4.1  VERIFICATION  OF  THEORY 

The  square  root  of  the  beam  degradation  of  the  RMS  value  of 
the  phase  aberation  as  given  earlier  is 

-  *j2  k  An  a  JUl  (11) 

Here  we  use  1  to  denote  the  RMS  value  of  the  phase  aberration;  k 
is  the  wave  number  equal  to  2*/\ .  Since 

t  =*  L/4  {12} 

as  measured  for  many  shear  and  wake  flows^)  we  get 


k  An 


a  L 


(  13) 


To  test  the  analytical  expression  above,  we  have  measured 
for  single  mixing  layers,  with  An  varying  from  0  to  2.45  x  10“4. 

To  avoid  large  velocity  nonuniformity  in  the  free  stream,  a  150 
mesh  screen  was  positioned  over  the  orifice  plate  {see  Section 
4.3).  Interf erograros  for  phase  aberration  studies  were  5.08  cm 
wide,  and  20  cm  high  in  the  stream-wise  direction.  Interf erogr ams 
for  mixing  layer  growth  studies  were  7.1  cm  wide  and  20  cm  high. 
Figures  8a  through  8f  are  typical  inte rf erogr ams  obtained. 


For  each  phase  aberration  interf erogr am ,  a  middle  portion  was 
digitized  and  processed  to  obtain  the  magnitude  of  phase  aberra¬ 
tion.  The  processed  area  was  7.2  cm  high  and  4.1  cm  wide  and  the 
lower  edge  of  the  area  was  7.2  cm  downstream  of  the  edge  of  the 
splitter  plate.  Along  the  length  of  this  area,  the  mixing  layer 
thickness  increased  by  about  50%.  The  average  thickness  of  layers 
studied  varied  from  0.6  to  0.9  cm.  The  experimental  results  of  14 
layers  studied  are  listed  in  Table  1.  Typical  phase  aberration 
contours  of  the  analyzed  inte rf eregrams  are  shown  in  Figures  9a 
and  9b.  These  show  the  magnitude  of  the  phase  aberration  normal¬ 
ized  to  2»  for  the  interferometer  laser  wavelength  of  5260  A. 

(All  other  phase  aberration  data  presented  in  this  section  are  for 
■  =  5360  &.)  The  field  of  view  of  the  inter ferograms  extends 
upstream  to  the  edge  of  the  splitter  plate,  0.63  cm  from  the  screen 
surface  (see  Fig.  3b). 


6.  fiatt,  k.G. ,  "Turbulent  Mixing  or  Passive  and  Chemically 
Reacting  Species  in  a  Low-Speed  Shear  Layer,"  J.F.M., 
Vol.  82,  1977,  pp.  53-95. 
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igure  Ga/b 


Beam  Degradation  (a),  Viewed  Perpendicular  to 
Mixing  Plane,  and  Mixing  Layer  Growth  (b) ,  Viewed 
Parallel  to  Mixing  P^ne,  for  In  =  6.2  x  10”5, 

PN  =  30  psig,  P  =  28  psig  and  \  =  5260  H 

(See  Table  1 ) 
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Beam 

for 


Degradation  (e)  and  Mixing  Layer  Growth  (f) 
An  =  24.5  x  10” P  =  30  psig,  P  =  7  psig 


-  5263  R  (See  Table  1) 
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TABLE 


Figure  9 


1  IHCH  *  0.403  W4Vt 
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/»***  ( b ) 


Phase  Aberration  Isometric  (a)  and  Iso-Contour  (b) 
Maps  for  An  =  18.6  x  10-5,  PN  =  30  psig,  PM  =  16 

psig  and  X  =  5260  £  (See  Table  1} 
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The  functional  dependence  of  <£'  on  An  can  be  seen  in  Figure 
10.  Here  the  values  of  phase  aberration  per  unit  mixing  layer 
thickness  ^’/2tL  are  given.  We  can  see  that  /  2tL  increases 
linearly  with  An.  The  solid  line  plots  the  theoretical  values  of 
*/2rL  computed  from  Eq.  (13)  with  a  =  18%.  This  value  of  a  is 
in  accordance  with  measured  values  obtained  in  a  large  number  of 
shear  and  wake  flows. (6)  Agreement  between  experimental  and 
analytical  results  can  be  seen. 

To  study  how  varied  with  the  streamwise  distance  x,  sev¬ 
eral  inte rf erograms  were  divided  into  five  contiguous  portions  and 
analyzed;  each  portion  was  a  rectangle  4  cm  high  and  4.5  cm  wide. 
The  results  are  listed  in  Table  2  and  plotted  in  Figure  11.  The 
values  of  L  and  x  were  measured  at  the  center  of  each  portion. 

The  results  were  obtained  for  two  values  of  An,  2.45  x  10~4  and 
1.86  x  10~4 .  (For  the  ABEI,  laser.  An  =  1.82  x  10-4) .  It  can 
be  seen  that  L  and  had  finite  values  at  x  =  0  and  then 
increased  linearly  with  x.  (See  the  following  subsection  for 
further  discussions  on  the  finite  values  of  L  at  x  =  0) .  The 
solid  lines  for  present  its  theoretical  values  based  on  the 
measured  values  of  L  and  with  a  -  14%.  This  value  of  a  is 
slightly  lower  than  that  obtained  from  Figure  10,  but  both  are  in 
agreement  with  measured  values  in  other  mixing  layer  experi¬ 
ments. 

In  summary,  the  data  obtained  in  this  series  of  experiments 
indicate  that  is  proportional  to  An  and  L  as  predicted  by 
Eq.  (2)  or  Eq.  (13)  . 

4.2  EFFECT  OF  VELOCITY  SHEAR 

The  growth  of  the  mixing  layer  thickness  depends  both  on  the 
intensity  and  macroscale  of  the  turbulence  in  the  free  stream  as 
well  as  the  velocity  difference  between  the  two  streams.  Brown 
and  Roshko^4)  also  reported  a  slight  dependence  of  the  growth 
rate  on  the  ratio  of  the  densities.  In  this  series  of  experi¬ 
ments,  we  studied  the  growth  rate  of  the  mixing  layer  between 
nitrogen  and  four  different  gas  mixtures.  A  150-mesh  screen  was 
positioned  over  the  orifice  plate  to  produce  uniform  free-stream 
velocity  (see  Section  4.3).  In  each  stream  the  turbulence  inten¬ 
sity  was  2%  at  1  cm  downstream  of  the  screen. 

Typical  interf erogram  results  are  shown  in  Figures  12a 
through  12c.  It  can  be  seen  that  L  grows  with  downstream 
distance  x  and  the  growth  rate  increases  with  higher  shear.  At 
x  =  0,  the  layer  thickness  L  is  about  1.5  mm  for  all  mixing  layers 
studied.  The  virtual  origin  of  the  layer  where  L  =  0  is  thus 
extrapolated  to  be  several  centimeters  upstream  of  the  edge  of  the 
splitter  plate.  The  layer  thickness  is  measured  with  the  method 
discussed  in  Section  3.1. 
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Evolution 

Thickness 

of  Phase  Aberration  and  Mixing  Layer 
Along  Stream-wise  Direction;  X  =  5260  R 
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Sabine(7)  and  Abramovich(8)  suggested  that  the  growth 

rate  of  a  shear  layer  thickness  is  a  function  of  the  ratio 


R 


+ 


u 

u 


2 

2 


(14) 


where  subscript  "1"  denotes  the  lighter  gas  and  "2"  the  heavier 
gas.  In  Figure  13  we  thus  plot  the  growth  rate, 

L'  =  dL/dx 


against  the  shear  ratio,  R. 

The  data,  which  cover  the  region  of  practical  interest,  |R|< 
0.4,  are  also  presented  in  Table  3.  L'  can  be  observed  to 
increase  with  the  absolute  value  of  R.  For  small  values  of 
[R |  (<  0.1)  where  shear  layers  in  aerowindows  are  likely  to  operate 
the  growth  rate  is  about  5%.  Measurements  were  taken  for  four 
density  ratios  of  1.33,  2.40,  3.18  and  7.0  (with  respective  values 
of  An  of  0.72  x  10“4,  1.66  x  10-4,  1.96  x  10"4  and  2.45  x 
10“4) ,  but  no  marked  effect  of  the  density  on  growth  rate  was 
observed. 

In  Figure  13  the  shear  flow  data  in  this  experiment  are  pre¬ 
sented  together  with  data  of  Brown  and  Roshko. (45  Their  data 
show  that  for  values  of  R  higher  than  0.4,  the  growth  rate 
increased  linearly  with  R;  a  more  marked  effect  of  density  on 
growth  rate  was  also  reported.  Our  data  in  the  region  of  interest 
of  | R|  <0.4  appear  to  converge  to  the  lowest  of  the  growth  rates 
reported  by  Brown  and  Roshko.  It  should  be  noted,  however,  that 
the  Brown  and  Roshko  data  were  obtained  by  a  somewhat  subjective 
visual  observation  of  the  shadowgraphs. 

The  Brown  and  Roshko  coherent  structures  could  also  be 
observed  in  our  experiments  (see  Figure  12c) .  These  structures 
become  much  more  pronounced  tor  1 R  j  >0.3.  It  is,  thus  advisable 
for  aerowindow  application  to  use  shear  layers  with  |r  |  < 0.1. 


7.  Sabin,  C.M.,  "An  Analytical  and  Experimental  Study  of  the 
Plane,  Incompressible,  Turbulent  Free-Shear  Layer  with 
Arbitrary  Velocity  Ratio  and  Pressure  Gradient,"  Trans.  ASME. 
D.  87,  1965,  pp.  421. 

8.  Abramovich,  G.N. ,  The  Theory  of  Turbulent  Jets,  M.I.T.  Press 
1963. 
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Figure  13  Effect  of  Velocity  Shear  on  Mixing  Layer  Growth 

Rate;  Lines  Through  Brown  and  Roshko  Data  Points 
(Ref.  4)  are  Sabin  and  Abramovich  (Refs.  7  and  8) 
Semiempirical  Theory  for  Large  1 R I 
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The  effect  of  shear  on  phase  aberration  is  shown  in  Figure 
14.  The  aberrations  were  obtained  for  the  rectangular  middle  por¬ 
tion  of  the  shear  layers  which  was  7.2  cm  high,  4.1  cm  wide  and 
with  its  lower  edge  7.2  cm  downstream  of  the  splitter  plate  edge. 
Here,  the  normalized  values  of  ^*/2»  are  plotted  against  R. 

For  |  R j  <0.1  the  normalized  pnase  aberration  is  about  0.18  x  lO'* 
for  X  =  5260  A. 

4.3  EFFECT  OF  FLOW  PLATE 

In  this  series  of  experiments,  we  studied  the  effect  of  a 
screen  positioned  over  the  orifice  plate  on  the  free-stream  veloc¬ 
ity  as  well  as  the  mixing  layer  growth  rate  and  phase  aberration. 

The  two  flow  plates  are  shown  in  Figures  3a  and  3b.  The 
screen  used  is  a  150-mesh  screen  with  34%  open  area.  The  junction 
between  the  screen  and  splitter  plate  must  be  carefully  made. 
Various  designs  for  the  junction  between  the  screens  and  splitter 
plate  extension  were  tried.  Figure  15  shows  the  nonrandom  non¬ 
uniformity  introduced  in  the  mixing  layer  by  one  unsuccessful 
design,  illustrated  in  Figure  16.  Itie  design  finally  decided  upon 
(Figure  17}  gave  very  satisfactory  results.  The  design  uses  two 
splitter  plate  extensions  with  a  fine  mesh  screen  clamped  between 
them.  Clamping  insures  good  alignment  between  the  two  extensions 
and  introduces  minimal  flow  disturbance  near  the  splitter  plate. 
The  wires  of  the  screen  are  oriented  45°  with  respect  to  the 
splitter  plate  to  avoid  formation  of  large-scale  nonuniformities. 
Flow  near  the  splitter  plate  was  measured  by  a  hot-wire  along  the 
extension  with  the  wire  positioned  1  mm  away  from  the  extension 
and  3  ram  down  stream  of  the  screen.  Flow  uniformity  along  the 
splitter  plate  extension  was  comparable  to  that  in  the  free 
stream.  This  design  was  thus  used  for  the  experiments  done  in 
this  test. 


The  free-stream  velocity  profiles  for  the  screened  and 
unscreened  flow  plate  are  shown  in  Figure  18.  A  screen  of 
150-mesh  was  used.  Quite  noticeable  differences  between  the  two 
profiles  can  be  observed.  Without  a  screen,  the  flow  is  very  non- 
uniform.  The  scales  of  nonuniformity  are  observed  to 
(3  or  4}  times  as  large  as  the  orifice  spacing.  This 
result  of  jet  coalescence  developed  in  orifice  plates 
area  ratios  of  less  than  about  60%.  (9) 


be  several 
i  s  the 
with  o  pe  n 


9.  Kuchemann,  D.  and  Sterne,  L.H.G.,  (1964)  ,  Progress  in 

Aeronautical  Sciences,  V5;  “The  Design  of  Low  Speed  Wind 
Tunnels,"  p.  12,  13,  28,  29,  30,  31. 
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Figure  15 


Non  Random,  Non-Uniformity  introduced  by  the 
Unsuccessful  Flow  Plate  Modification  Illustrated 
in  Figure  16;  •  =  5260  R 
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Figure  16  Modified  Flow  Plate  of  Unsuccessful  Design 
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Figure  17  Modified  Flow  Plate  of  Successful  Design 
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J9496  TRANSVERSE  DISTANCE  (cm) 


18  Effect  of  Screening  on  Flow  Uniformity;  Measure¬ 
ments  Made  1.27  cm  Downstream  of  Orifice  Flow  Plate 
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To  obtain  uniform  flow  from  the  flow  plate,  it  is  important 
that  a  fine-mesh  screen  be  positioned  over  the  plate  so  that  jets 
from  the  plate  can  be  "broken"  up.  "Jet  Breaking"  is  more  pro¬ 
nounced  with  finer  mesh  screens  which  give  higher  flow  resist¬ 
ance.  Figure  19  shows  the  velocity  profiles  measured  for  various 
screens.  Screens  with  mesh  number  higher  than  100  give  good 
uniformity.  These  finer  screens  are  usually  manufactured  with  an 
open  area  ratio  <40%,  and  jet  coalescence  occurs.  The  scales  of 
nonuniformity,  as  observed  in  Figures  18  and  19,  are,  however, 
kept  small  due  to  the  small  mesh  size,  A  mesh  size  of  150  was 
chosen  to  be  used  in  our  experiment,  as  a  compromise  between 
strength  and  effectiveness  in  smoothing  out  free-stream  velocity. 

The  effect  of  the  screen  on  mixing  layer  growth  rate  and 
phase  aberration  was  studied  with  two  values  of  An  of  0.62  x 
10-4  and  1.03  x  10-4.  The  shear  ratio,  R,  of  the  mixing 
layers  was  taken  to  be  about  0.06.  It  should  be  noted  that,  in 
order  to  keep  farfield  beam  degradation  (due  to  shear  layer 
turbulence)  negligible,  } R |  should  be  less  than  0.1.  As  shown  in 
Figures  20a  and  20b,*  tne  mixing  layer  growth  rate  is  much  higher 
for  the  unscreened  flow  plate.  The  growth  rate,  as  well  as  the 
phase  aberration  for  the  middle  portion  of  the  mixing  layer  (7.2  x 
4.1  cm2  area  with  the  lower  edge  7.2  cm  from  the  splitter  plate 
end)  are  listed  in  Table  4  and  presented  in  Figure  21;  the  dotted 
lines  on  the  plot  are  theoretical  extrapolations  based  on  Eq . 

(2).  For  the  unscreened  plate,  the  growth  rate  and  aberration  are 
about  twice  as  large  as  those  for  the  screened  flow  plate.  This 
means  that  the  beam  degradation  for  the  screened  flow  plate  is 
only  about  one-fourth  of  the  degradation  for  the  unscreened  flow 
plate.  Screening  is  thus  a  very  good  method  for  reducing  beam 
degradation. 

Another  advantage  of  screening  is  that  the  phenomena  of 
mixing  layer  bending  and  deflection  observed  for  the  unscreened 
flow  plate,  does  not  occur.  This  phenomena  as  reported 
earlier occurs  for  the  unscreened  flow  plate  when  there  is  an 
appreciable  difference  between  the  pressures  of  the  two  plenum 
chambers  upstream  of  the  plate.  For  the  screened  plate,  similar 
pressure  differences  in  the  plenums  do  not  give  rise  to  notice¬ 
able  mixing  layer  bending.  (See,  for  example,  the  mixing  layer  in 
Figure  8f  where  the  pressures  in  the  two  plenums  are  7  psig  and  30 
psig.)  For  the  unscreened  flow  plate  the  mixing  layer  bending 
toward  the  gas  stream  with  high  plenum  pressure  was  most  severe 
within  a  few  centimeters  downstream  of  the  flow  plate  and  became 
unnoticeable  further  downstream. (2)  The  apparent  reason  for  the 
bending  is  that  the  gas  stream  with  higher  plenum  pressure  has  a 
higher  kinetic  energy  as  shown  in  Appendix  B;  the  gas  stream  with 
higher  kinetic  energy  will  create  a  negative  static  pressure  differ 
ential  and  "pull"  the  mixing  layer  toward  its  side.  This  "ejector 
pump  effect"  is  most  severe  when  the  difference  between 


*The  field  of  view  in  Figure  20a,  for  the  screened  flow  plate, 
extends  upstream  to  the  edge  of  the  splitter  plate,  0.63  cm  from 
the  screen;  the  field  of  view  in  Figure  20b,  for  the  unscreened 
flow  plate,  extends  upstream  to  the  flow  plate. 
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Effect  of  Mesh  Size  on  Flow  Uniformity;  Measure¬ 
ments  Made  1.27  cm  Downstream  of  Orifice  Flow  Plate 
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Figure  20  Effect  of  Screening  on  Mixing  Layer  Growth  for 

An  =  6.2  x  10~~>  and  =  5260  S  (see  Table  4)  with 
150  Mesh/Inch  Screen  (a)  and  Without  Screen  (b) 
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the  kinetic  energies,  and  thus  the  static  pressures,  of  the  two 
streams  is  highest.  For  the  unscreened  flow  plate  this  is  the 
region  near  the  flow  plate  where  jets  from  the  orifices  have  not 
yet  merged  to  form  a  uniform  flow.  As  the  orifice  plate  has  an 
open  area  ratio  of  40.8%,  the  kinetic  energies  of  the  gases,  and 
consequently  their  difference,  in  the  jetting  region  is  six  times 
higher  than  those  in  the  uniform  flow.  The  modified  flow  plate 
with  a  fine-mesh  screen  over  the  orifice  plate  has  effectively 
eliminated  this  jetting  region  above  the  screen  and  thus  reduced 
mixing  layer  bending  to  unnoticeable  levels. 

4.4  EFFECT  OF  EXPOSING  ONE  FREE  STREAM  TO  A  STAGNANT  REGION 

As  shown  in  Figure  1,  the  end  flow  in  the  ABEL  laser  is 
exposed  to  a  stagnant  region  which  is  connected  to  the  aero- 
window.  To  study  the  effect  of  this  stagnant  region  on  the  mixing 
layer,  experiments  were  performed  in  the  flow  setup  shown  in 
Figure  2  with  nitrogen  and  helium  flowing  from  the  two  plenum 
chambers.  Experiments  were  performed  under  unexposed  and  exposed 
conditions.  Under  the  unexposed  condition  both  streams  were  con¬ 
fined  by  the  chimney  walls.  Under  the  exposed  condition  the 
nitrogen  stream  was  exposed  to  the  stagnant  region  and  the  helium 
stream  confined  by  the  chimney  wall;  the  stagnant  region  was 
5.3  cm  away  from  the  mixing  layer,  and  its  height  he  was  varied 
from  10  cm  to  18  cm. 

Interferometer  results  are  shown  in  Fiqures  22  and  23.  In¬ 
spection  reveals  that  the  mixing  layers  in  the  exposed  condition 
were  not  markedly  different  than  those  under  unexposed  con¬ 
ditions.  No  instability,  such  as  wavering  of  the  layers  was 
observed.  Quantitative  comparisons  of  the  phase  aberration  are 
presented  in  Table  5.  The  results  for  layers  under  both  con¬ 
ditions  are  roughly  equal  within  experimental  error. 
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Figure  22  Mixing  Layer  Growth  with  One  Stream  Open  to  a 

Stagnant  Region  for  in  =  24.5  x  10-^,  he  =  10  cm, 

P  =  30  psig,  PM  =  7  psig  and  =  5260  R  (Compare 


is  due 


with  Figure  S f )  .  The  large-scale  distortion  is  due 
to  poor  optical  quality  Plexiglas  in  the  interferometer 
beam  path  (other  interf erograms  were  taken  through  a 
high  optical  quality  quartz  window) 
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TABLE  5 

PHASE  ABERRATION  *'/2ir  FOR  DIFFERENT  VALUES  OF  h 
FOR  MIXING  LAYERS  BETWEEN  He  AND  N2  FLOWING  AT  35.5eM/S 
AND  28.5  M/S  RESPECTIVELY* 


♦Values  of  aberration  were  obtained  from  analyzing  rectangular 
areas  4.5  cm  wide  and  7.2  cm  high;  x  is  the  streamwise  distance 
of  the  center  of  each  area;  \  -  5260  A 
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5.0  SUMMARY  AND  CONCLUSIONS 


The  following  conclusions  can  be  obtained  from  the  results  of 
this  present  test  series: 

1)  The  analytical  expression  for  random  beam  degradation 
with  laser  propagation  normal  to  two-dimensional 
turbulent  mixing  layers, 


B 


„  (  2  *  \2  2.2 
2  \~rj  “  ' 

is  verified  by  measurements. 

2)  From  the  expression  above,  we  see  that  three  methods  can 
be  adopted  to  reduce  beam  degradation: 

a)  reducing  L£  by  fluid  mechanical  means, 

b)  reducing  An  by  using  other  laser  gases  (see  Ref.  3, 
Appendix  E) , 

c)  Introducing  multiple  layers  (see  Ref.  3). 

3)  If  An  is  fixed,  fine-mesh  screens  positioned  over  the 
flow  plate  are  quite  effective  in  reducing  Li  .  The 
growth  of  the  mixing  layer  is  reduced  roughly  by  one-half 
with  the  use  of  a  150-mesh  screen.  This  means  beam 
degradation  is  reduced  by  a  factor  of  about  four.  Mixing 
layer  bending  is  also  avoided  by  using  the  screened  flow 
plate . 

4)  Higher  velocity  shear  increases  the  mixing  layer  growth. 
The  growth  rate  is  a  function  of 


The  growth  rate  for  a  layer  with  |r|  ^0.1  is  about  5%  if 
a  screened  flow  plate  is  used. 

5}  Exposure  of  a  free  stream  to  a  stagnant  region 

sufficiently  far  away  from  the  mixing  layer  does  not  have 
an  appreciable  effect  on  the  mixing  layer  characteristics. 

Beam  degradation  due  to  turbulent  mixing  for  the  ABEL  device 
with  An  -  1.82  x  10-4  can  now  be  calculated.  The  results  for  a 
screened  flow  plate  as  in  Figure  11,  with  An  =  1.86  x  10"4  and 
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R  =  0>G6,  can  be  used  directly.  For  the  same  value  of  R,  if  the 
mixing  layer  height  in  the  flow  direction  is  30  cm,  the  average 
beam  degradation*  is  0.133  x  (2*)^  for  x  =  5260  A.  In  the  case 
of  the  ABEL  COj  laser  with  X  =  10.6  ^m,  the  beam  degradation  for 
each  mixing  layer  is 


£1 

I 


(0.133) 


x 


( 2  *■) 


x 


/ 0.526 
\  10.6 


1.3% 


The  total  degradation  of  the  laser  beam  due  to  the  shear  layers  is 
thus  about  4.1%  since  the  beam  passes  1.6  times  through  the  two 
end  layers. (3)  This  random  degradation  appears  sufficiently 
small  to  be  acceptable. 


*For  shear  layers,  the  average  beam  degradation  between  xj.  and  X2 
is  =  ('J’j  +  i>22  +  ^  A )  /  3  where  subscripts  "I"  and  "2"  denote 

values  at  and  respectively. 
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APPENDIX  A 


THE  OPTICAL  SCIENCES  CENTER  CODE 
FOR  INTERFEROGRAM  ANALYSIS,  "FRINGE" 


The  FRINGE  code  is  an  inte rf erogram  analysis  program  capable 
of  handling  either  straight  or  circular  fringes.  Straight  fringes 
do  not  have  multiple  values  when  a  line  is  drawn  nearly  perpendi¬ 
cular  to  the  fringe.  The  circular  fringes  have  no  multiple  values 
along  a  radius  from  the  pupil  center  or  center  of  symmetry. 

The  program  reads  in  an  array  of  (x,y)  points  representing 
the  fringe  maximum  or  minimum  and  has  mnemonic  codes  to  tell  the 
computer  when  fringes  are  missing  (i.e.,  on  edge  of  pupil  or  in 
central  obstruction) .  The  program  has  a  very  flexible  and  easily 
used  input  control  format  and  output  option  selection.  Internally 
the  program  takes  the  (x,y)  array  points  and  first  fits  them  to  a 
polynomial  of  the  form 

f(x,y)  =  Ax  +  By  +  C  +  D(x2  +  y2)  (A-l) 

The  first  three  terms  represent  the  tilt  while  the  last  term 
give  the  defocus.  The  magnitude  of  D  represents  the  defocus  in 
waves  at  the  edge  of  the  pupil.  The  fit  (x,y)  is  obtained  via  a 
least-squares  routine  and  represents  the  fringe  order.  The  dif¬ 
ference  between  the  fit  and  the  stored  array  of  data  yields  the 
phase  errors,  i.e., 


OPDi  =  ORDi  -  f  (x  , y)  (A-2) 

where  OPDi  is  the  optical  path  difference  at  point  i  in  the 
matrix  and  ORDi  is  the  fringe  order  at  point  i  and  f(xi,  y ^ ) 
is  the  value  of  the  polynomial  fit  to  the  data. 

For  straight  fringes,  aberration  terms  such  as  astigmatism 
and  third-order  coma  can  be  added  to  the  polynomial  either  sep¬ 
arately  or  together.  If  this  option  is  selected,  the  phase  errors 
or  residuals  from  the  least-squares  fit  will  be  relative  to  this 
higher-order  polynomial.  For  the  circular  fringe  case,  tilt  and 
defocus  are  removed  and,  if  desired,  polynomial  terms  up  to  10th 
order  are  fit  to  the  data.  The  conic  constants  can  be  solved  for 
or  specified  as  desired.  The  program  has  an  option  to  average 
several  interf erograms  or  subject  a  baseline.  This  is  important 
if  the  effects  of  the  measuring  devices  are  to  be  removed. 
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There  are  several  plotting  options.  Currently  implemented 
are  three-dimensional  isometric  plots  and  contour  maps.  Line 
printer  contour  maps  are  also  obtainable  but  not  used. 

In  summary,  this  program  has  the  capability  to  analyze  inter- 
ferograms.  Phase  error  plots  are  produced  and  polynomial  coeffi¬ 
cients  representing  wave  aberrations  are  generated.  Both  straight 
and  circular  fringes  can  be  handled. 
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APPENDIX  B 


RELATION  BETWEEN  KINETIC  ENERGY  AND  PLENUM  PRESSURE 

The  kinetic  energy  of  a  flowing  stream  in  the  laser  cavity  is 

PV2  =  pa 2  M2  (B- 1 ) 

where 

pi  gas  density  in  the  cavity 

a:  sonic  velocity  of  the  gas  in  the  cavity 

M:  Mach  number  of  the  flowing  gas  in  the  cavity. 

But  from  perfect  gas  law,  we  have: 

pH2  =  )  (>RT)  =  y  p  (B-2 ) 

where 

p:  cavity  pressure 
R:  gas  constant 
T:  cavity  temperature 
Y:  ratio  of  specific  heat  capacities 

From  Eqs.  (1)  and  (2),  we  get 


(B-3 ) 


Now,  the  flow  in  the  sonic  orifice  flow  plate  is  given  by 


(B-4) 
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where 


V:  gas  velocity  in  the  cavity 

A:  total  flow  plate  area  for  uniform  flow  far  from  the  flow 

plate  or  effective  jetting  area  for  jetting  flow  near  the 
unscreened  flow  plate. 

A*:  sonic  orifice  area 


7  +1) 
2(>-l) 


K:  defined  as  y'Y+l 

PQ :  plenum  pressure 
T0 :  gas  temperature  in  the  plenum. 

For  low  Mach  Number,  we  have 

T = T0  so  that  from  Eg.  (4)  we  get 


M 


From  Eqs.  (3)  and  (5)  we  get 


pV 


2 


(B-5) 


(B-6) 


which  is  our  relation  between  kinetic  energy  and  plenum  pressure. 
We  see  that  the  kinetic  energy,  at  low  Mach  number  ,  is  approx¬ 
imately  independent  of  the  gas  temperature.  For  the  unscreened 
ABEL  flow  plate  the  effective  jetting  area  is  40.8%  of  the  total 
flow  plate  area,  so  the  kinetic  energy  in  the  jetting  region  is 
about  six  times  higher  than  the  kinetic  energy  in  the  uniform  flow 
region. 

The  analysis  above  applied  for  choked  flow  across  the  orifice 
plate  and  the  condition  for  choked  flow  is 
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